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ABSTRACT. The behaviour of resonances in the spin-orbit coupling in Celestial Me- 
chanics is investigated. We introduce a Hamiltonian nearly-integrable model describing 
an approximation of the spin-orbit interaction. A parametric representation of periodic 
orbits is presented. We provide explicit formulae to compute the Taylor series expansion 
in the perturbing parameter of the function describing this parametrization. Then we 
compute approximately the radius of convergence providing an indication of the stability 
of the periodic orbit. This quantity is used to describe the different probabilities of capture 
into resonance. In particular, we notice that for low values of the orbital eccentricity the 
only significative resonance is the synchronous one. Higher order resonances (including 
1:2, 3:2, 2:1) appear only as the orbital eccentricity is increased. 
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§1. INTRODUCTION 



The dynamics of coupling between revolutional and rotational motions of a satellite is 
investigated. In particular, we consider an oblate satellite moving on a Keplerian orbit 
around a central planet and rotating about an internal spin-axis (sec Celletti 1990, Celletti 
1994, Celletti and Falcolini 1992, Celletti and Chierchia 1998, Goldreich and Peale 1966, 
1970, Henrard 1985, Murdock 1978, Peale 1973, Wisdom 1987 and references therein for 
related papers on this subject). Whenever the ratio between the revolutional and rotational 
periods is rational, say = - for some positive integers », o, a spin-orbit resonance 

rot Q 

(of order p : q) is encountered. Astronomical observations show that most of the evolved 
satellites or planets of the solar system are trapped in a 1:1 or "synchronous" spin-orbit 
resonance. The most familiar example is provided by the Moon; due to the equality of the 
periods of rotation and revolution, the Moon always points the same face to the Earth. The 
only known exception is provided by Mercury, which moves in a 3:2 resonance, namely after 
two revolutions around the Sun, Mercury makes three rotations about its spin-axis. An 
intriguing question is the explanation of the different occurrence between plenty resonances 
(like the 1:1 or 3:2) and empty resonances (corresponding to any other different ratio of 
the periods). In this paper we try to address this problem as follows. We introduce a 
mathematical non-autonomous, one-dimensional model describing a physically relevant 
approximation of the spin-orbit problem (see §2). The equation of motion takes the form 

X - e fx{x,t) = , 

for a suitable function / depending also on the orbital eccentricity e of the Keplerian orbit; 
the parameter e is proportional to the equatorial oblateness coefficient of the satellite. 

In §3 we introduce a parametric representation of a periodic orbit with frequency | as 

x = 9 + u{9) 

y = v{d) deT = R/(27rZ) , 

where u, v are suitable functions depending analytically on e and 9{t) — 9{0) + ^t. We 
provide explicit formulae for the computation of the coefficients of the Taylor series ex- 
pansion around £ = of the function u, say u{9) = YljLi '^ji^)^'' • Finally, we compute an 
approximation of the radius of convergence 

£p(-) = lim Y- 

(for a fixed 6 = ^o)- This quantity provides a measure of the regularity of the periodic 
orbit and is used to provide a qualitative argument toward a higher capture probability 

of the synchronous resonance. We remark that a recent method based on a quantitative 
version of the Implicit Function Theorem has been implemented in (Celletti and Chierchia, 
1998) to construct Birkhoff periodic orbits. 

The results and conclusions are presented in §4. The graphs of £p(|) vs. | show that for 
the Moon and Mercury, the 1:1 and 3:2 resonances respectively seem to be the most suitable 
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ending-states. Moreover we explore the birth of resonances as the orbital eccentricity is 

increased. In particular, for low values of e the synchronous resonance is the only one 
present. As the eccentricity increases, higher order resonances do appear. We relate these 
results to the evolutionary history of the satellites. 

We remark that an exhaustive explanation of this scenario requires the consideration of 
dissipative forces and their role in stabilizing the resonances. We demand this problem to 
a future study. 



§2. THE SPIN-ORBIT MODEL 

We briefly introduce a mathematical model describing the "spin-orbit" interaction in Ce- 
lestial Mechanics as follows. 

Let S he a triaxial ellipsoidal satellite orbiting around a central planet P. Let Trev and 
Trot be the periods of revolution of the satellite around P and the period of rotation about 
an internal spin-axis. A p : q spin-orbit resonance occurs whenever 



J- rot Q 

In particular, when p — q — 1 we speak of 1:1 or synchronous spin-orbit resonance, which 
implies that the satellite always points the same face to the host planet. In the solar 
system most of the evolved satellites or planets (like, e.g., the Moon) are trapped in a 
1:1 resonance (Astronomical Almanac, 1997). The only exception is provided by Mercury 
which is observed in a nearly 3:2 resonance. 

We introduce a mathematical model describing the spin-orbit interaction. In partic- 
ular we assume that 

i) the center of mass of the satellite moves on a Kcplcrian orbit aroinid P with semimajor 
axis a and eccentricity e (i.e., we neglect secular perturbations on the orbital parameters); 
a) the spin-axis is perpendicular to the orbit plane (i.e., we neglect the "obliquity"); 
in) the spin-axis coincides with the shortest physical axis (i.e., the axis whose moment of 
inertia is largest); 

iv) dissipative effects as well as perturbations due to other planets or satellites are ne- 
glected. 

Let A < B < C he the principal moments of inertia of the satellite. We denote by r and 
/, respectively, the instantaneous orbital radius and the true anomaly of the Kcplcrian 
orbit. Finally let x he the angle between the longest axis of the ellipsoid and the periapsis 
line (see Figure 1). The equation of motion governing the spin-orbit model under the 
assumptions i) — iv) may be derived from the standard Euler's equations for rigid body. 
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In normalized units (i.e. assuming that the mean motion is one, 2TT/Trev = 1) we obtain: 

X + e{-f sin(2a;-2/) = 0, (2.1) 
r 

where e = | is proportional to the equatorial oblateness coefficient (and the dot 

denotes time diff'crentiation) . Notice that (2.1) is trivially integrated when A = B or m 
the case of zero orbital eccentricity. 

A ^^p : q periodic orbit" (or "Birkhoff periodic orbit with rotation number p/q") is a 
solution of (2.1) such that 

x{t + 2'Kq) = x(t) + 2ttp ; 

the above relation implies that after q orbital revolutions around the central body the 
satellite makes p rotations about the spin-axis. 



Due to assumption i), the quantities r and / are known Keplerian functions of the time; 
therefore we can expand (2.1) in Fourier series as 



oo 

X + e ^ ^(?,e) sui{2x-mt) = 0, (2.2) 

m^O,m= — oo 

where the coefficients W{Y,e) decay as powers of the orbital eccentricity as Vr(^,e) oc 
el"T^-2| (^ggg Cayley 1859, for explicit expressions). 

A further simplification of the model is performed as follows. According to assumption 
iv), we neglected dissipative forces and any gravitational attraction beside that of the 
central planet. The most important contribution comes from the non-rigidity of the satel- 
lite, which provokes a tidal torque due to internal friction. Since the magnitude of the 
dissipative effects is small compared to the gravitational term, we simplify (2.2) retaining 
only those terms which are of the same order or bigger than the average effect of the tidal 
torque. Finally we consider an equation of the form 

N2 

X + e ^ ^(^^e) sin(2a; — mt) = , 

m^O,m=Ni 

where A^i and A^2 are suitable integers (depending on the structural and orbital properties 
of the satellite), while W^(y,e) are truncations of the coefficients Vr(^,e) (which are 
power series in e). 

According to the above discussion, in the case of the Moon-Earth system we consider the 
following equation of motion: 

^ + £ [(-| + Y6^ sm{2x-t) + 

+ (1 - + ^e^) sin(2a; - 2t) + ( - ^e^) sin(2a; - 3t) + 

2 lb 2 lb (2 3") 

,17 2 115 4, . , , ,845 o 32525 . , , ^ ' ^ 

+ (— e^) sin(2a; - At) + ( e^) sin(2a; - 5t) + 

T V 2 Q ' ^ ^ ^ 48 768 ' ^ ' 



+ ^e^sin(2x-bt) + ^l^e^ sin(2x - 7t) ] = 0, 
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having taken A^i = 1 and A^2 = 7 in (2.2). In what follows we shall always assume that 
the leading equation is provided by (2.3), though in principle we should consider different 
truncation orders (A'^i, A^2) depending on the intrinsic parameters of the satellite. For 
example, in the Mercury-Sun case the above criterion would suggest to take A^i = —17 
and N2 = 6. However we checked that our results do not change significantly as new 
terms are added to eq. (2.3). Precisely, more terms alter the outputs for any initial 
condition, provoking a global rescaling of the results, without conditioning the qualitative 
conclusions we shall draw about the behaviour of the resonances. As an illustration we 
report in Figure 2 the Poincare map around the synchronous resonance for e = 0.004, 
e = 0.2. The resonance is surrounded by librational surfaces, whose amplitude increases as 
the chaotic separatrix is approached. Outside this region rotational surfaces can be found. 



§3. PARAMETRIC REPRESENTATION OF PERIODIC ORBITS 

In this section we provide a parametric representation of the solution associated to a 
periodic orbit with frequency u; = |. 
We rewrite eq. (2.3) in compact form as 

X - efx{x,t) = , 

namely 

X = y 

A periodic orbit with frequency a; = | is defined parametrically by the set of equations 

x = e + u{9) 

y^l + Du{e) eeT , ^^-^^ 

where u — u{9) is a suitable continuous function, depending analytically on £, with the 
property that the flow in the ^-coordinate is linear, i.e. ^ — > ^ + |t after a time t; the 
operator D acts on a function u = u{9) as 



Inserting (3.2) in the system (3.1) one obtains the second order differential equation 

L>V^)-£/,(^ + «(^),-) = 0. (3.3) 
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Due to the analyticity in the perturbing parameter s of the function u, we can expand u 
in Taylor series around £ = as 



oo 

«(^) = ^«„(^)£" , (3.4) 



for suitable p-periodic functions Un{0). Let us expand the terms Un{0) in Fourier series as 



k 

where the coefficients u'^^ (and therefore Un{0)) can be constructed explicitely as follows. 
Write the function fx{x,t) as 

Ux,t) = I fm (e^(2— t)_e-i(2.-mt)) ^ ^35) 

where in the case of eq. (2.3) the set M. consists of the integers 1,2,...,7 and the fm can 
be easily identified by comparison with eq. (2.3). From eq. (3.5) it follows that 

k=±l,m&M 

Define a power series for the exponential as 

oo 
n=0 

It can be easily verified that the coefficients 6^"^^ satisfy the following recursive relations: 

h=l 

Inserting (3.4) in (3.3) one gets 

oo oo 

J2D'M0)e'^ = iY'"" E ^M-U- ((3-6) 

n=l n=l k=±l,m&M 

Moreover, since an explicit computation of the generic term in the summation at the l.h.s. 
of (3.6) provides 

^ k 



a comparison of the above relation with the r.h.s. of (3.6) yields the explicit expression of 
the coefficients u^^^ ■ 

Recasting the above formulae, the function u{9) can be written as 

CO „ 

u{e) = 5^«n(^)£" , With un{e) = 5^4"^sin(nfe-) , 

n=l k ^ 



for suitable Fourier indexes n/-; the fact that UniO) is a sum of sines is due to the specific 
form of eq. (2.3) which contains only sine-terms. 

As a measure of the regularity of the function tt(^), we compute for any d = 9q the radius 
of convergence of the Taylor series as 

e {I) = lim (3.7) 

(which in fact seems not to depend on the choice of ^o)- The radius of convergence provides 
an indication of the stability of the periodic orbits, showing an approximate value of the 
perturbing parameter e — ep{^) at which the transition from elliptic to hyperbolic periodic 
orbits takes place. A numerical investigation of the phase space suggests that as e grows 
the periodic orbit is surrounded by librational curves with increasing amplitude, until they 
leave place to a chaotic regime as e approaches £p(|)- We make use of this property in 
order to study the behaviour of the periodic orbits as the frequency | is varied. 



§4. RESULTS AND CONCLUSIONS 

According to standard evolutionary theories, satellites and planets were rotating fast in 
the past; a constant decrease of the period of rotation about the spin-axis was provoked by 
tidal friction due to the internal non-rigidity. Therefore, a common scenario suggests that 
celestial bodies were slowed down until they reached their actual dynamical configuration, 
being typically trapped in a 1:1 or 3:2 (for Mercury alone) resonance. This hypothesis im- 
plies that higher order resonances (i.e., 2:1, 5:4, 7:3, etc.) were bypassed during the slowing 
process. However there is actually no convincing explanation concerning the mechanism 
of escape or capture into a spin-orbit resonance. In this section we argue that there is a 
greater probability of capture into the synchronous or 3:2 resonance. More precisely, we 
compute an approximate value of the radius of convergence of the parametric representa- 
tion of periodic orbits which provides, as remarked in §3, a measure of their stability. We 
let the eccentricity vary in a reasonable (astronomical) range of values. In particular, we 
consider the following set of periodic orbits with frequencies | where 

9=1, ...,14, p = 9+ 1,...,15 , 
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and 

p=l,...,U, q = p+l,...,15 . 

In order to have a better precision around the main resonances we consider also the periodic 
orbits with frequencies 

1 1 1 

± —-^ , 1± 
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where k = 1, 10. We computed the quantity — ^ for j = 1,...,30, providing a 

good approximation to the radius of convergence. The choice of the point 9 — Oq does not 
influence the result, since the radius of convergence seems to be generally independent on 
the initial condition. 

As an example we report in Figure 3 the radius of convergence £p{j) as a function of the 
number of iterations j = 1, 30 (in this case was set equal to 1.56). 

For a fixed value of the eccentricity we compute £p{^) corresponding to the set of rational 
numbers ^ listed above. Figures 4 and 5 show the graphs of £p(^) versus | for, respec- 
tively, the Moon's eccentricity (i.e., e = 0.0549) and Mercury's eccentricity (e = 0.2056). 
In order to have a more detailed inspection of the behaviour around the resonances of 
astronomical interest, we report in panel a) the results around the 1:2 resonance, while 
panel b) corresponds to the synchronous resonance, panel c) to the 3:2 resonance and panel 
d) to the 2:1 resonance. 

A comparison between Figures 4 and 5 suggests that locally the 1:1 and 3:2 resonances have 
a higher probability of capture. In fact, the periodic orbits around these commensurabilities 
have very low radii of convergence indicating that the most stable orbits correspond to 
the periods ^ = i or £ = |. Moreover, comparing Figure 46, 4c with Figure 56, 5c 

q L ^ q z 

we notice that there is a slightly preference for the Moon to end-up in the synchronous 
resonance, while Mercury might be trapped in the 1:1 or 3:2 resonance, in agreement with 
the astronomical observations. 

This discussion leads to the question of the existence of attracting tori for a dissipative 
system (i.e., including tidal forces) corresponding to the most stable resonances. We plan 
to address this problem in a later study. 



Next we analyze the behaviour of the stability of periodic orbits for different values of the 
eccentricity. Figure 6 shows £p{^) versus ^ for e = 0.001 (Figure 6a), e = 0.01 (Figure 66), 
e = 0.06 (Figure 6c), e = 0.2 (Figure 6d). For low values of the eccentricity (Figure 6a) 
there is only one marked resonance corresponding to the synchronous commensurability. 
Moreover, the amplitude of the curve around the 1:1 resonance provides a measure of 
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the size of the region of hbrational motion surrounding the resonance. For values of | 

bigger than |, the quantity ep(|) increases indefinitely with |. As e = 0.01 (Figure 66), 
beside the 1:1 resonance there appears the 3:2, while the 1:2 and 2:1 resonances are already 
present, but with very small amplitudes of librational motion, which become meaningful 
as e = 0.06 (Figure 6c). For this value of the eccentricity we have the complete set of main 
resonances, i.e. 1:2, 1:1, 3:2, 2:1. Minor resonances appear as e is increased up to e = 0.2 
(Figure 6d), where the 5:4 and 7:4 resonances become evident. 

These results indicate the existence of a strict relation between the value of the eccentric- 
ity and the birth of resonances. Satellites with low eccentricity are suitable candidates for 
ending-up in the synchronous resonance. However, even when different resonances arise, 
the 1:1 periodic orbit seems the most likely final state due to the amplitude of librational 
regime around it. When the eccentricity is highly increased, several new resonances ap- 
pear. In particular, for Mercury's eccentricity (compare with Figure 6d) different fates are 
possible, including the 3:2 resonance. 

As a final remark, we suggest that a further development of this study would include the 
effect of dissipative terms and their role in driving satellites to select their ending states. 
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FIGURE CAPTIONS 



Figure 1: The spin-orbit geometry. 

Figure 2: Poincare map associated to eq. (2.3) around the synchronous resonance in 
the (a;, i;)-plane for e = 0.004, e = 0.2. The initial conditions are {x,y) = (0,0.8), 
(0,1.5), (0,1), (0,1.1), (1.8,1), (2.1,1), (1.57,1). Notice that eq. (2.3) is 7r-periodic in 
X. Librational and rotational regions are divided by the chaotic separatrix. 

Figure 3: The radius of convergence Sp{j) vs. the degree of approximation of the limit 
(3.7), j = 1, 30. Here = 1-56 and e = 0.0549. 

Figure 4: The radius of convergence £p(|) vs. | for e = 0.0549 (i.e., the eccentricity of 
the Moon). The dots correspond to the actual computations associated to the frequencies 
listed in the text. Lines are due to interpolation performed by the graphic program. In 
abscissa are reported the rotation numbers ^ (the 1:1 resonance corresponds to 1, the 1:2 
to 0.5 and so on). 

a) 1:2 resonance, b) 1:1 resonance, c) 3:2 resonance, d) 2:1 resonance. 

Figure 5: The radius of convergence £p(|) vs. | for e = 0.2056 (i.e., the eccentricity of 
Mercury). The dots correspond to the actual computations associated to the frequencies 
listed in the text. Lines are due to interpolation performed by the graphic program. In 
abscissa are reported the rotation numbers | (the 1:1 resonance corresponds to 1, the 1:2 
to 0.5 and so on). 

a) 1:2 resonance, b) 1:1 resonance, c) 3:2 resonance, d) 2:1 resonance. 

Figure 6: Plot of the radius of convergence £p(^) vs. | for the frequencies listed in the 
text, a) e = 0.001, b) e = 0.01, c) e = 0.06, d) e = 0.2. 
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